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Electronic structure and transport properties of the fullerene §2 and the metallofullerene Gd@§2 
are investigated with density functional theory and the Landauer-Buttiker formalism. The ground 
state structure of Gd@§2 is found to have the Gd atom below the C-C bond on the C2 molecular axis 
of §2. Insertion of Gd into §2 deforms the carbon chain in the vicinity of the Gd atoms. Significant 
overlap of the electron distribution is found between Gd and the §2 cage, with the transferred Gd 
electron density localized mainly on the nearest carbon atoms. This charge localization reduces 
some of the conducting channels for the transport, causing a reduction in the conductivity of the 
Gd@§2 species relative to the empty §2 molecule. The electron transport across the metallofullerene 
is found to be insensitive to the spin state of the Gd atom. 

I. INTRODUCTION 

Endohedral metallofuUerencesi have attracted wide interest due to their functional characteristics and poten- 
tial applications in the field of nanomaterials and biomedical sciencei^ Recently, metallofullerene doped nanotubes 
("peapods") have also attracted experimental attention due to their structural and electronics properties j^i^ and 
have been proposed as a possible self-assembled quantum computing architecture J- Additionally, recent experimental 
. studies of Gd@§2 in single wall carbon nanotube ((Gd@§2)@SWNT) peapods have shown novel transport behaviori 
5t . The first step in understanding the properties of the peapod structures is an adequate treatment of the individual 
Gd@§2 fuUerenes. 

In earlier Gd@§2 theoretical calculations by Kobayashi and Nagase^^ the relaxation of the fullerene cage was not 
considered, nor the possibility of novel geometries for the Gd atom in the fullerene cage was not considered; the pre- 
dicted ground state spin multiplicity M = 9, resulting from ferromagnetic coupling between the Gd f-electrons and the 
odd electron on the fullerene cage, is in contrast to electron spin resonance determinations of a M=7 ground state, due 
to antiferromagnetic coupling of — 1.8cm~^ii2i The M=7 ground state is also supported by magnetic studies piSiiLiSii^ 
In order to resolve this discrepancy, we present results of density functional calculations for the equilibrium geometries, 
magnetic and transport properties of Gd@§2, made with the Landauer-Buttiker transport formalism. Section Hll gives 
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a brief description of the theoretical procedure, Section HTll present our results, and Section Hvl draws conclusions and 
discusses possible future directions. 

II. COMPUTATIONAL METHODS 

The equilibrium geometry and the total energy of Gd@§2 are calculated using density functional theory (DFT)ii4 
One of the primary considerations involved in these calculations is determination of a suitable basis set and exchange 
correlation functional. We therefore compared the bond length, binding energy and spin multiplicity values of Gd 
and Gd dimer (Gd2) optimized using several different basis sets and exchange correlation functionals, to experimental 
valueglS and theoretical all-electron relativistic calculations''^^. As we discuss in Section Fill Al the results are quite 
sensitive to the choices of both basis set and exchange correlation functional. A relativistic effective core potential 
CEP-121G^^ for the atoms with the generalized gradient approximation of Becke's exchange functionaU^ and Lee- 
Yang-Parr correlational functional^^ (B3LYP) gave results comparable to experiment, and were thus used for the 
subsequent calculations. 

The total energies and forces are calculated using the linear combination of atomic orbitals (LCAO) molecular orbital 
approachfi^ We employed a 6-31G basis for the carbon atoms and a CEP-121G basis set for Gd. The computations 
were performed using GaussianOS for various spin multiplicitiesn^S The geometries were optimized without symmetry 
constraints by minimizing the total energy and requiring the forces to vanish within a 10""^ a. u. /Bohr threshold, at 
every atom site. To find the ground state structure for Gd@§2, initial positions of the Gd atom inside the §2 cage were 
generated, then optimized without constraining any of the atomic coordinates. The calculations were repeated for 
different spin multiplicities to determine the ground state spin configuration, as we discuss in Section [ill Bl Although 
the partitioning of charge to specific atoms or groups within a molecule is not uniquely defined within the postulates 
of quantum mechanicSfSi a recent experimental evaluation of density functional charge schemes?^ found the commonly 
used MuUikcn charge analysis^"'^ to be generally deficient for charge calculations, and Natural Population Analysis^^ 
(NPA) to be the most generally accurate method of those studied. 

To study the transport, we created a model comprised of §2 located between two Au contacts, with and without 
the endohedral Gd atom. The numbers of Au atoms on each side of the contact were varied to study the transport 
dependence on contact geometry. Fig.l shows one of these structures with a triangular contact of 3 Au layers in 
which the layer adjacent to §2 has one Au atom, the second has 3 Au atoms and the third layer has 6 Au atoms. 
For calculation of electron transport through Gd@§2 and §2, we used the non-equilibrium Green's function based 
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Landauer-Buttiker formalisini2^i2^iS& Neglecting spin-flip processes, the total current due to coherent scattering is 
given by Ispin-coherent^I"' + , whcre /" and are the contribution to current from spin up (a) and down (/3) 
states, respectively. The spin specific contributions are expressed as the integral over the injection energy of the 
tunnelling electron, E, 

taken between the electrochemical potentials fix and fj,2, 

Ml = Ef -{q,V/iqi+q2)) (2) 
M2 = Ef + {q2V/{qi+q2)), (3) 

where qi and 52 are charge accumulated in the left and right contacts respective, V is the applied bias, and Ef is the 
Fermi energy of the gold contact, taken as -5.53 eV for our calculation. f{E,fi) is the Fermi distribution function and 
T"^^^(E, V) is the transmission function that represents the sum of transmission probabilities for electrons of a given 
spin through the contact/metallofuUerene complex, obtained from the overlap and Hamiltonian matrices determined 
by the DFT calculation. The explicity included Au atoms in our model system described above are used to obtain the 
coupling matrices for calculation of self-energy functions. We have used the local density of states of the 6s-band 
of bulk gold (0.035 eV per electron spin) to approximate the Green's function of the (bulk) Au contacti2& 

III. RESULTS 
A. Properties of Gd and Gd2 

First, we discuss the geometry, ionization potential, and electron affinity of Gd and Gd2 and compare these with 
experiment. This comparison is provided to assess the accuracy of our theoretical procedure. 

To assess how well our method accounts for the properties of the Gd atom, we have calculated the ionization potential 
of Gd and Gd2 as well as the binding energy and Gd2 dimer bond length. The ground state Gd atom was found to 
have a spin multiplicity of 9 for all functionals that we have considered with the CEP-121G effective core potential 
basis functions.il Resuhs for the LSDA^ BSLYPi^ and EWQlS"^ functionals are presented in Supplemental Table 
I. The calculated binding energy value (1.7595 eV) and bond-length (2.926 A) for Gd2 with CEP-121G as frozen core 
basis function and B3LYP as functional for spin multiplicity of 19 agrees well with the experimental binding energy 
fl.784±0.35 eVU^ and the relativistic all-electron theoretical bond-lenrth f2.895 The ionization Dotential for 
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Gd and Gd2 were calculated as 5.014 eV and 4.119 eV respectively. Based on these results, the CEP-121G basis set 
for Gd and 6-31G basis set for C atoms, together with the B3LYP functional, were used for the Gd@C82 calculations. 

Although we have employed a frozen core, the Gd 5s and 5p core electrons are allowed to relax in the molecular 
calculations. Atomic calculations show that the relaxation of these states is often crucial for an accurate determination 
of valence state properties because of the overlap of the 5d electron orbitals with the other n=5 orbitals. 

B. Structure of Gd@§2 

We have taken a §2 fuUerene cage with point group symmetry ^2{G2v) and optimized the structure with Gd 
inside at various positions, as described in Sectional In the §2(02,;) fuUerene cage, the C2-axis goes through the 
center of one six-membered ring and one C-C double bond. Our calculations reveal that the ground state structure 
of Gd@§2 has the Gd atom situated adjacent to the C-C double bond on the symmetry axis, as shown in Fig. 2. 
This is in contrast to reported results for other metal atoms such as Sc, Y and La in ^2fSi22iS2i^2i2ii^ where the 
metal atoms are found to be centered inside the fuUerene cage on top of the C2v axis six-membered luig^^ Recent 
synchrotron radiation powder structure analysis experiments of Gd@§2 are consistent with our calculated placement 
of the Gd atom.^^ The distances between Gd and the two nearest C atoms are 2.38 A and 2.41 A respectively. The 
binding energy of Gd to §2 is found to be 5.6435 eV (130.1 kcal/mole), comparable to those for the Sc, Y, and La 
metallofuUerenes;* the energy levels of §2 and Gd@§2 are shown in Supplemental Figure 1. Addition of Gd into §2 
slightly lengthens the proximal C-C bond (from 1.425 A to 1.473 A) atom, as well as the neighboring C-C bonds (by 
0.01 to 0.04 A). This deformation results in an energy difference between the §2(C2y) symmetry structure and the 
Gd@§2 optimized deformed structure when recalculated without Gd of 0.109 eV (2.51 kcal/mole). 

As discussed in Section ^ the NPA charge partitioning tends to give better results; our Gd@§2 calculations agree 
with this analysis, as the MuUiken charge analysis indicates charge transfer of 1.43 electrons from Gd atom to the 
C82 cage, and NPA indicates a charge transfer of 2.43 electrons. The latter value agrees well with electron energy 
loss spectroscopy (EELS) experiments.'^^ The electrons are localized near the two proximal carbon atoms, resulting 
in strong charge density overlap between the Gd atom and these C atoms. NPA analysis also gives a Gd d-orbital 
population of 0.48 electrons, consistent with a Dewar-Chatt back-bonding interaction between the Gd atom and the 
carbon double bondu^^ The deformation of the C-C bonds near the Gd is also consistent with back-bonding, as the 
back donation of charge from the C-C bonds to the d-orbitals weakens (and hence lengthens) the bonds. Additionally, 
the Gd /-orbital population of 7.02 is consistent with the charge transfer from the cage of 0.04 electrons empirically 



5 

determined by Nadai et al. for reproducing X-ray magnetic circular dichroism experimentsi^ 

In order to study the difference in Gd placement, as compared to the Sc, Y, and La cases, we performed single point 
calculations, in which the Gd atom was placed in the La-like position, adjusting the distance between the hexagon and 
the Gd atom to account for the smaller radius of the Gd^'^ ion as compared to the La^'^ ion;^ The lowest total energy 
structure we found occured when the Gd atom was placed 2.10 A above the plane of the hexagon; this structure had 
a total energy of 0.104 hartrees (2.83 eV) higher than the ground position determined by our optimizations. NPA 
analysis determined nearly identical charge distribution on the Gd atom in both cases, attesting to the importance 
of Gd-C back-bonding in stabilizing the structure. However, an elementary electrostatic analysis, using the NPA 
atomic charges as point ions, gives a Gd-cage potential energy of -1.015 and -0.9366 hartrees for the ground and 
La-like positions respectively; the difference between these nearly accounts for the difference in total energies. This 
interpretation is consistent with resonant photoemission spectroscopy experiments which suggest that the Gd-C bond 
is primarily ionic in character^ 

We have also examined the effect of different spin configurations on the optimized the Gd@§2 structure. The 
ground state structure was found to have 7 unpaired electrons on Gd and one unpaired electron on the §2 cage 
aligned antiparallel to the Gd spins, with total spin multiplicity M=:7, in agreement with magnetic studiesiiSiiiiiSii^ 
The next lowest energy structure has M=9, with seven unpaired electrons at Gd and the odd §2 cage electron 
oriented parallel to the Gd spins. The energy difference between these two configurations is 2.6 mcV (16.1 cm^^), 
which compares favorably to the experimentally determined value of 14.4 cm~^^ The detailed energies with various 
spin multiplicities are listed in Table I. 

C. Electronic transport 

In order to see the effects of contact geometry, we have taken between one and three layers of Au atoms to represent 
the contact. To obtain the contacts, single Au atoms were placed on each side of the Gd@§2 molecule, shown 
schematically in Fig. 1, at a distance of 2.015 A from the carbon atoms, along the C2v symmetry axis. To build the 
larger contacts, the bond lengths and geometry of bulk Au were used to add 3 Au atoms to the second layer and 6 Au 
atoms to the third layer, as shown in Fig. 1. Although it is known that the conductance can be strongly dependent on 
the contact structure, ■^^ this transport model serves as the basis for understanding the qualitative transport properties 
of Gd@§2. 

Calculated I-V characteristics for §2 and Gd@§2 with various spin states and Au contact geometries are shown in 
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Figs. 3, 4, 5, and Supplemental Fig. 3. The ground state (M=7) transport with a single Au atom contact on each side 
of the fuUerene is shown in Fig. 3. Addition of the Gd atom into the §2 cage leads to a reduction in the current. To 
understand this difference, we computed the density of state (DOS) for both structures, and found the metallofullerene 
to have a reduced DOS near the Fermi energy as compared to §2. These results are shown in Supplemental Fig. 2. 
The I-V characteristics for the two- and three-layer contact structures were found to be qualitatively the same as for 
the single Au atom contact case, and are shown in Supplemental Fig. 3 and Fig. 4, respectively. Quantitatively, 
the current is seen to increase with the addition of further gold layers, due to the increasing number of conduction 
channels. 

Calculations for the M=9 structure, which is structurally degenerate to the M=7 system, are shown in Fig. 5. We 
find only a small difference in current between the two spin configurations, and the conduction properties of spin up 
and spin down electrons (as calculated in Eq.JQ)) were found to be insensitive to the Gd spin configurations, suggesting 
that the spin of the Gd is shielded by the cage carbon atoms and that transport occurs consequently only through 
the fuUerene cage. This is also supported by plots of the highest occcupied and lowest unoccupied molecular orbitals 
(HOMO/LUMO), shown in Figs. 6 and 7, which indicate that conduction occurs primarily through the §2 cage and 
shows no signature of conducting paths involving the Gd atom. This suggests that incorporation of Gd decreases the 
number of available conduction channels between the contacts, as a result of the charge localization effect noted in 
Section ImBl 

By attaching the Au contacts, one expects some charge transfer from the metallic contact to the fuUerene at zero 
bias, but this effect (which is the same for both structures) is relatively small due to the weak coupling between the 
fuUerene and metal contact. Using the NPA method described in the previous section, the charge transfer from the 
left and right Au contacts to the §2 cage were found to be 0.09 and 0.11 electrons for single layer of Au, 0.12 and 
0.13 electrons for the two layer contact and 0.13 and 0.15 electrons for the three layer contact. Use of these charge 
values in the expression for the electrochemical potential, Eq.Q, had no effect on the calculated transport properties. 
Despite the fact that Gd is closer to one of the electrodes, we observed no diode like behavior; we attribute this to 
the the conduction occuring primarily through the fuUerene cage, which is approximately symmetrical with respect 
to both contacts. 
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IV. CONCLUSIONS 

We have analyzed the energetics, structure, and transport properties of Gd@§2. The present resuhs indicate that 
the CEP-121G basis set for Gd, with the B3LYP density functional for exchange and correlation, is a satisfactory 
combination for DFT calculations for large Gd endofuUerene systems. Our calculated structure for Gd@§2 confirms 
the geometry determined by x-ray diffraction data, indicating that the Gd is located adjacent to the C-C double 
bond on the symmetry axis, in contrast to the other Group-3 metallofuUerenes. Furthermore this geometry gives 
calculated spin configurations that agree with ESR and magnetic experiments, in contrast to previous theoretical 
studies. Using the Landauer-Buttiker formalism for transport, we find that conduction occurs primarily through 
the §2 cage, and that charge donation from the Gd atom to the cage disrupts these conduction channels. We find 
no evidence for electron-spin dependent transport effects due to the spin state of the Gd atom. Besides providing 
a solid methodological and structural basis for future calculations of (Gd@§2)@SWNT structures, the absence of 
HOMO /LUMO density on the Gd atom — and thus relative lack of interaction with the surrounding nanotube — helps 
explain why the M4 and M5 peak edges of the (Gd@§2)@SWNT EELS spectrum were found to be identical to that 
for Gd@§2, indicating that the Gd valence state is unaffected by the surrounding nanotubeii 
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TABLE I: Calculated total energies and binding energies for the optimized structures of §2 and Gd@§2 as a function of Gd 
spin multiplicity, M. The M=7 state of Gd@§2 is found to be 2.6 meV below the M=9 state in total energy. 



Atom{M) 


TE{au) 


BE{eV) 


§2(C2„) 


-3123.805927 




Gd(a§2(M=5) 


-3235.8783289 


4.4066 


Gd@§2(M=7) 


-3235.9237819 


5.6435 


Gd@§2(M=9) 


-3235.9236904 


5.6409 


Gd@§2(M=ll) 


-3235.87890696 


4.4224 




Au Gd@C82 



FIG. 1: Gold contact model used for transport calculations across Gd@§2 (see part b). The first layer adjacent to §2 has one 
Au atom, the second has 3 Au atoms and the third layer has 6 Au atoms. 
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FIG. 2: Optimized structure of Gd@§2. The ground state structure of Gd@§2 has the Gd atom (purple) situated adjacent to 
the C-C double bond on the C2v symmetry axis that connects a six-membered ring with a C=C bond on the other side of the 
fuUerene. The distances between Gd and the two nearest C atoms are 2.38 A and 2.41 A respectively. The lines between Gd 
and C are used to indicate the closest carbon atoms. 
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Voltage (V) 

FIG. 3; I-V plot for single Au atom contacts. Total conduction (sum of spin up and spin down conduction) for Gd@§2 is lower 
at all bias voltages than for pure fuUerene §2. Conduction of spin up (alpha) and spin down (beta) electrons is similar at all 
bias values. 
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FIG. 4: I-V plot for contact geometry shown in Figure 1. Note the qualitative similarity to the other contact geometries (Figs. 
3, Supplement Fig. 2). 
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FIG. 5: I-V plot for conduction in the spin multiplicity M=7 ground state, and the M=9 statee with energy 2.6 meV higher in 
energy. Conduction characteristics for spin polarized and unpolarized currents are similar for both cases, with neither showing 
strong spin-dependent transport effects. This suggests that the current is mainly carried through the §2 cage, and does not 
involve the Gd atom. 
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FIG. 6: Lowest Unoccupied Molecular Orbital (LUMO) plot of Gd@§2 showing localization on the §2 cage without significant 
LUMO orbitals on the Gd atom. Blue and green are used to indicate the positive and negative sign of the wavefunction, 
respectively. The purple sphere indicated the Gd atom. 




FIG. 7: Highest Occupied Molecular Orbital (HOMO) plot of Gd@§2. Blue and green are used to indicate the positive and 
negative sign of the wavefunction. HOMO orbitals are localized on cage of §2 and there is no HOMO orbital on Gd (purple 
sphere). From both HOMO and LUMO plots, it is evident that the current is mainly carried through the §2 cage, as opposed 
to the Gd atom. 



